home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Info-Mac 4
/
Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso
/
Science
/
MAL ©P.f.Howden 1⁄1⁄89
/
INTEGRATE
< prev
next >
Wrap
Text File
|
1989-07-12
|
1KB
|
17 lines
1 CLS:PRINT"INTEGRATOR=Area (quadrature) or Line (arc) for single-value Functions.":INPUT "MENU: 0=A Function:1=Area under N DATA points:=";M:DEFDBL A-I,O-Z:IF M=0 THEN 4
2 CLEAR:DEFDBL A-I,O-Z:INPUT "ORDINATE SPACING H=";H:INPUT "N(ODD)=";N:DIM X(N):FOR J=1 TO N:PRINT "INPUT X";J;:INPUT "=";X(J):F=F+X(J)*(3+(-1)^J):NEXT J
3 PRINT "AREA INTEGRAL=";(F-X(1)-X(N))*H/3:PRINT"END :":PRINT:GOTO 2
4 PRINT "TYPE FUNCTION ON LINE 10. eg: 10 F=X^3-2":INPUT "MENU: 0=AREA under curve: 1=LINE integral: ?";L
5 IF L=1 THEN INPUT "DIFFERENTIAL (.0001?)=";Z
6 GOTO 20
7 IF L=1 THEN X=X+Z:GOSUB 10:F1=F:X=X-2*Z:GOSUB 10:F2=F:X=X+Z:F=SQR(1+((F1-F2)/(2*Z))^2):RETURN
10 F=3*X^5-4*X^3
19 RETURN
20 INPUT "LOW LIMIT=";A:INPUT "TOP LIMIT=";B:INPUT "ACCURACY (.001?)=";D
30 H=(B-A)/2:X=A:GOSUB 7:P=F:X=B:GOSUB 7:P=H*(P+F):I1=3*P:N=1
40 S=0:FOR K=1 TO N:X=A+(2*K-1)*H:GOSUB 7:S=S+F:NEXT K
50 I0=P+4*H*S:IF ABS(I0-I1)<D THEN 70
60 P=(P+I0)/4:N=2*N:H=H/2:I1=I0:GOTO 40
70 PRINT "INTEGRAL=";I0/3:PRINT"END :":PRINT:GOTO 5
MACINTOSH LISTING